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ABSTRACT 

We present XMM-Newton observations of the Chandra -detected nuclear X- 
ray source in NGC 4561. The hard X-ray spectrum can be described by a model 
composed of an absorbed power-law with T = 2.51^3, anc ^ column density N H = 
1.91q2 x 10 22 atoms cm -2 . The absorption corrected luminosity of the source is 
L(0.2 - 10.0 keV) = 2.5 x 10 41 ergs s" 1 , with bolometric luminosity over 3 x 10 42 
ergs s _1 . Based on the spectrum and the luminosity, we identify the nuclear X- 
ray source in NGC 4561 to be an AGN, with a black hole of mass Mbh > 2 x 10 4 
M Q . The presence of a supermassive black hole at the center of this bulge- 
less galaxy shows that black hole masses are not necessarily related to bulge 
properties, contrary to the general belief. Observations such as these call into 
question several theoretical models of BH-galaxy co-evolution that are based 
on merger-driven BH growth; secular processes clearly play an important role. 
Several emission lines are detected in the soft X-ray spectrum of the source which 
can be well parametrized by an absorbed diffuse thermal plasma with non-solar 
abundances of some heavy elements. Similar soft X-ray emission is observed in 
spectra of Seyfert 2 galaxies and low luminosity AGNs, suggesting an origin in 
the circumnuclear plasma. 

Subject headings: galaxies: active — galaxies: individual (NGC 4561) — galaxies: 
nuclei — X-rays: individual (NGC 4561) 
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INTRODUCTION 



The past decade has seen extraordinary growth in our understanding of supermassive 
black holes (S MBHs), with secure det ections, mass measurements and new demographic in- 
formation (see iFerrarese fc Fordll2005l and references therein). Knowledge of the mass func- 
tion of SMBHs directly affects our understanding of SMBH formation and growth, nuclear 
activity, and the relation of SM BHs to the format ion and evolution of galaxies in hierarchical 
cold dark matter models (e.g., iMenci et al.l 12004 ). The cumulative mass function needed to 
explain the energetics of hi gh redshift quasars implies that all galax ies in the local universe 



should host a SMBH (e.g.. Marconi et all 12004 IShankar et al.ll2004f ). 



Observationally, however, we do not know whether every galaxy hosts a SMBH. Tradi- 
tional methods of finding SMBHs, viz. stellar dynamics and gas dynamics are powerful only 
at the high-mass end of the SMBH mass function: the BH sphere of influence cannot be 
resolved for BH masses less than 10 6 M beyond a distance of a couple of Mpc, even with 



HST (IFerrarese Fordll2005l ). One way, and perhaps the most efficient way to find SMBHs 
in galaxies is to look for active SMBHs. In fact, looking for AGN activity is perhaps the 
only viable way to probe the low-mass end of the local SMBH mass function. X-ray observa- 
tions provide the best opportunity for this purpose because X-ray emission is an ubiquitous 
property of AGNs and X-rays can penetrate obscuring material that might hide an AGN 
at other wavelengths. Indeed, X-ray observations have detected AGN activity in what were 
thought to be "no rmal" galaxies in clusters (e.g., iMartini et al.ll2002l ) and in the field (e.g., 
Brand et al.ll2005h . 



In an effort to study the demographics of local SMBHs, we had undertaken a Chan- 
dra program to look for nuclear X-ray sources in nearby optically "normal" spiral galaxies 
(within 20Mpc). This program was highly su ccessful; we discoy e red AGNs in the nuclei of 



what were th ought to be "normal" galaxies (IGhosh et al.l 120081 . iGhoshl 120091 . iGhosh et al. 
2011 ; see also Zhang et al. 2009 ; Desroches fc Ho! 2009 ). The nuclear X-ray sources, how 



ever, could be stars, binaries, supernova remnants or AGNs. Through extensive spectral, 
timing, and multiwavelength analysis we classified the nuclear X-ray sources and found 17 
(out of 56 surveyed) that are almost certainly low-luminosity AGNs. Thus at least 30% of 
"normal" galaxies are actually active. The inferred luminosities of these sources range from 
10 375 to 10 42 ergs s _1 . In a few objects where SMBH masses were known from stellar/gas 
velocity dispersion methods, we find accretion rates as low as 10~ 5 of t he Eddington limit, 
comparable to what has been found in LINERS (e.g JPudik et al.ll2005| ). 



We then expanded upon the initial Chandra survey by usi ng the sample of SINGS 
galaxies (Spitzer INfrared Galaxy Survey; iKennicutt et al.l 120031 ). Out of the 75 SINGS 
galaxies, 60 have data in the Chandra archive and we detected nuclear X-ray sources in 36 
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of them (IGrier et al.l 120111 ). This once again shows that X-ray observations are far more 
efficient at detecting AGNs than optical. 

As noted above, through multiwavelength analysis we have shown that a large fraction of 
the Chandra -detected nuclear X-ray sources are indeed AGNs. Additional ly, using statistica l 
arguments we have shown most of them to be AGNs (Ghosh et al. 2010. lGrier et al.l 1201 lh . 
We obtained XMM-Newton spectra of several of our Chandra -detected nuclear X-ray sources 
with the goal of obtaining secure identifications as either AGNs or other contaminants. In 
this paper we focus on the bulge-less galaxy NGC 4561 for the following reason. 

The mass of the supermassive black holes (BHs) in centers of galaxies was found to be 



corre lated with t he bulge luminosity o f host galaxies (M B h Bulge relation; iMagorrian et al. 



19981 ; revised in iGiiltekin et al.ll2009l ). Even a tighter correlation was late r found between 



the BH mass and the velocity dispersion (a) of the bulge ( M B h - ° relation; iGebhardt et al. 



2000l . iFerrarese fc Merrittl |2000| . iMerritt fc Ferraresd 120011 ) ■ Basically, the mas s of the black 
hole seems to be correlated with the mass of the bulge (IHaring fc Rixi l2004f) . The above 



relations for norm al galaxies also extend to active galaxies (e.g., iMcLure fc Dunlopl 12002 



Woo fc Urryj|2002l ). These results were interpreted to imply that the formation and growth of 
the nuclear black hole and the bulge in a galaxy are intimately related, and several theoretical 
models have attempted to explain the ob served Mbh - & and Mbh -^Buige relations (e.g. 



Adams et al. 2001. Pi Matteo et al 



such as those of iHopkins et al. 



2006 



20031 ). The hydrodynamic cosmological simulations, 



naturally account for BH-galaxy co-evolution. In all 
such models, the bulge determines the nuclear BH mass. In the models of Volonteri and 
Natarajan (2009) seed BH masses are correlated to the host dark matter properties- thes e 
seeds could have been formed by direct collapse of pre-galactic halos (IBegelman et al.ll2006l ). 
Major mergers then trigger simultaneous BH growth and star formation resulting in a tight 
coupling between the two. A clear prediction of these models is that low-mass bulge-less 
galaxies today are unlikely to host nuclear BHs (Volonteri, Natarajan & G:ultekin 2011). 
Finding AGNs in bulge-less galaxies would certainly be a challenge for such models. Perhaps 
the BHs in bulge-less galaxies represent the seed BHs that have not yet grown. It has also 
been shown that the merger and SMBH growing process m ay not create a bulge b ecause of 
star formation and, mainly, because of supernova feedback (IGovernato et al.ll2010l ). Recent 



studies have shown that even mode rate-luminosity AGN s up to z~3 are powered mostly 



through internally driven proce sses (IMullanev et al. 



major role in triggering AGNs (ICisternas et al 



2011 



201 i f), and that merger s do not play a 



Schawinski et al.ll201ll ). Thus the role 



of mergers in the growth of SMBHs remains a matter of debate and SMBHs in bulge-less 
galaxies provide an important piece of the puzzle. Enlarging the sample of these objects 
is thus crucial to learn their properties and to be able to understand their formation and 
evolution mechanisms. 
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In this paper we present XMM-Newton data of the nuclear X-ray source in the bulge-less 
galaxy NGC 4561 which shows evidence of being an AGN. This secure identification of a 
SMBH in a bulge-less galaxy shows that a bulge is not necessary for the existence of a BH 
and the presence of SMBHs in bulge-less galaxies is more common than expected. 



2. NGC 4561 

NGC 4561 (Sdm) is a late-type bulge-less spiral galaxy at z=0. 00469. It follows the 
selection criteria used in our previous Chandra survey (close to face on with inclination less 
than 35° to ensure that the nuclear source is not obscured by the disk of the host galaxy; 
Galactic latitude \b\ > 30° to avoid obscuration and contamination from our own gala xy; and 



no kn own starburst or AGN activity) in which a nuclear X-r ay source was detected (|Ghosh 
2009 ). The optical spectra of NGC 4561 were analyzed by Kirhakos fc Steiner 199ol . who 



classified it as an H II region-like galaxy. 

The nuclear X-ray source was detected in a Chandra observation with 103 net counts 
with a count rate of 0.029 ± 0.003 ct/s. The source was found to be hard (HR = — 0. 53^13)- 
The spectrum was fitted with a power law T = 1.5 ± 0.3 and no intrinsic absorption (N# < 
1.7 x 10 21 cm" 2 ). With this model the flux was F(0.3 - 8 keV) = 2.5 x 10~ 13 ergs cm" 2 s" 1 
which corresponds to a luminosity of L(0.3 - 8 keV) ~ 5 x 10 39 ergs s -1 . This luminosity is 
a lower limit: the hardness of the source indicates that the source is likely to be absorbed. 
The quality of the Chandra data, however, was not good enough to determine the spectral 
shape accurately. The hardnes s of the emis sion and the high enough luminosity make this 



source a good candidate AGN (IGhoshl 120091 ) . 



The Chandra observations also showed a second source at 7" from the nucleus (source 
"B" here onward), which is soft (HR = —0.9) and its flux is only 10% of the nuclear source 
flux. Even though source B won't be resolved by XMM-Newton , it's effect on the spectrum 
of the nuclear X-ray source should be minimal. 



3. Observations and Data Reduction 

Our target was observed with XMM-Newton on the 2009 July 10. For the European 
Photon Imaging Camera (EPIC), the exposure times were 57022 s for MOS1, 57038 s for 
MOS2 and 55960 s for pn; all of these exposures were obtained using the thin filter in 
extended full frame. 

The data were processed and filtered using SAS v9.0.0 using tasks epchain and emchain. 
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Before the source extraction, we applied standard and temporal filters to the event list. 
For the standard filter we selected the energy to be between 0.2 and 15 keV for pn and 
between 0.2 and 12 keV for MOS, single or double pixel events (PATTERN <= 4) for 
pn and single, double or triple pixel events for MOS (PATTERN <= 12), with good flag 
values (#XMMEA_EP) for pn and (#XMMEA_EM) for MOS. For the temporal filters we 
created light curves for the 10-12 keV band selecting only events with PATTERN==0 (as 
recommended at the XMM SAS User Guide)Q. The good time intervals were created by 
rejecting the intervals with count rates higher than 0.5 cts/s for pn and 0.25 cts/s for MOS. 
The effective exposure time is 26.3 ks (for pn) corresponding to 47% of the observation time. 
The effective exposure times for MOS after applying the GTI are 36.7ks (64%) for MOS1 
and 37.6ks (66%) for MOS2. 

Within the circular extraction region of 20" radius the total number of counts detected 
in the PN data was 6638 cts with a rate of 0.282 ± 0.004 cts s _1 . The exposure time at the 
source (after vignetting correction) was 21.91 ks. Our source was detected in soft and hard 
X-rays with 4529.4 cts (0.2-2.5 keV) and 1581.7 cts (2.5-10.0 keV), which corresponds to 
HR = —0.48 ±0.01. The background level was of 2.13 cts/pixel. There were no considerable 
variations on the count rate during the observation. 

3.1. The spectra 

For extracting the spectrum we used a selection area of a 20" radius circle centered 
at the source and selected only events with FLAG==0 in order to obtain a good quality 
spectrum. When extracting the background spectrum we used a square of 41" side. Finally 
the RMF and ARF files were created and backscale was ran. Corrections for Out of time 
(OoT) events or Pile up were not needed. 

The spectrum was analyzed using Xspec vl2.6.0 (HEASOFT). We binned the spectra 
using grppha with a minimum of 50 cts per bin for the PN spectrum and 40 cts per bin for 
the MOS1 and MOS2 spectra. 

3.1.1. Checking for possible contamination 

We checked if spectrum was contaminated. This was a possibility given that there is a 
third source at 32" from the nucleus (source C here onward), which could contaminate the 



http: / /xmm.esa.int / external / xmm_user .support /documentation/sas_usg/USG / 
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source spectrum. The source extraction radius of 26.65" encircles 90% of the energy, so some 
contamination from source C is expected. In order to see how important this contamination 
was, we tried different selection areas when extracting the spectrum: (1) a circle with a 
radius of 20" centered at the source; (2) a circle with a radius of 26.65" (which encircles 90% 
of the total energy); (3) the same 26.65" circle but excluding a 26.65" circle around source 
C; and (4) a 26.65" radius semicircle on the opposite side of the source C. 

These selection shown in figure [TJ Emission line-like features were present in 

all of the 4 extracted spectra, discussed further in §4. Also, all of the spectra have similar 
values for T and for the observed flux showing that the contamination is negligible, and that 
the emission features present in the spectrum are real. 

3.1.2. The fitting process 

The spectral fitting was performed on the PN spectrum for energies E>0.3 keV. Because 
of the considerably lower signal to noise of the MOS data, it did not help to constrain the 
parameters any better. Once the final models were chosen, we tried them on the MOS1 and 
MOS2 spectrum and confirmed that they are consistent. 

To fit the spectra, different models were tried, always with a fixed galactic absorption 
of Nh = 2.11 x 10 20 cm -2 and a free intrinsic absorption. Both absorption components were 
represented with the "wabs" model for photoelectric absorption. 

First we tried with a power law, which by itself does not provide a good fit (x 2 =1.28 
for 114 dof) and overestimates the counts for high energies. To avoid the latter issue we 
fitted the absorbed power law using only the hard part of the spectra (E> 2.5 keV). To 
get a good fit in the hard range of energies, intrinsic absorption is needed (as suspected 
from the Chandra observation). We needed another component to model the soft-band, 
since the power law is almost completely absorbed in the soft-band (see figure [2]). The 
resulting photon index was T=2.46. From here on, when fitting the spectrum over the entire 
energy range, the T parameter was kept frozen at this value. In order to characterize the 
unresolved nearby source (source B), we also added a black body component with a flux of 
F(0.3-8keV)=2.5xlO -14 ergs cm~ 2 s _1 corresponding to the flux observed with Chandra and 
a temperature of kT=0.2 keV. Both the temperature and the normalization were always 
fixed. 

Since the hard-band is well fit by a power-law, our next task is to fit the soft component 
with different models. Because the spectrum shows "emission- like" features, we tried an 
absorbed "mekal" model, which characterizes emission from hot diffuse gas, but the fit was 
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not good (x 2 =2.35 for 112 dof, even worse than the absorbed power law with A% 2 = —117.3, 
figure |3]). This fit was done fixing the abundance to the Solar value. Leaving this parameter 
free does not improve the fit considerably (x 2 ,=2.27 for 111 dof, figure H]). 

We also tried to fit the soft component with the "vmekal" model, same as "mekal" but 
allows the abundances of each element to vary individually. A good fit was finally found using 
vmekal when leaving the abundances of some of the heavier elements as free parameters. This 
model has x 2 =l.ll for 107 dof (figure [SJ, that corresponds to A% 2 = 144.43 when compared 
to the fit with solar abundances. 

It is possible that in the soft-band we are observing the continuum reflected off some 
nearby scattering material. For this reason, we also tried the "reflionx" model (reflection 
by a constant density illuminated atmosphere). This model does reproduce some emission 
lines, but the fit is not good and it is also a worse fit than the absorbed power law (x 2 =1.39 
for 113 dof, Ax 2 = —11.18 when compared to the absorbed power law model, figure [6]). 



4. Results 

The best fitted model is composed of an absorbed power law with T = 2.5±8;| and N H = 
2.01q2 x cm ~ 2 an d an absorbed thermal plasma component with non-solar abundances 
for some of the heavy elements (O = 0.33^, Na = 34^\ Si < 0.2, Fe = 0.12+°;^, Ni 
< 0.4) with a temperature of kT = 0.59±g;g| keV and N H = 7t\ x 10 20 cm" 2 (and a black 
body for the source B). The fit to the spectrum is presented in figure |5]and the corresponding 
theoretical model in figure [3 The observed flux is F(0.2 - 10 keV) = 1.2 x 10~ 12 ergs cm -2 
s _1 which corresponds to a luminosity of L(0.2 - 10.0 keV) = 5.8 x 10 40 ergs s" 1 . The 
unabsorbed luminosity is L(0.2 - 10.0 keV) = 2.5 x 10 41 ergs s -1 . The bolometric luminosity 
of the source is therefore about 3.5 x 10 42 ergs s" 1 (a ssuming a bolometric correction factor 



of 14 for this luminosity; IVasudevan fe Fabianl 120071 ). putting it squarely in the luminosity 
range observed for AGNs. 



5. Discussion 



The hard-band spectrum of the nuclear X-ray source in NGC 4561 can be described as 
a power law with photon index V = 2.51q; 3 . Most of the soft component of the power law is 
absorbed (intrinsic absorption of Nh = 2.0loJ x 10 22 cm~ 2 ). The source shows a hardness 
ratio of HR = —0.48 ± 0.01 and an absorption corrected luminosity of L(0.2 - 10.0 keV) 
= 2.5 x 10 41 ergs s . The photon index and the high luminosity of the source are indications 
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of the presence of an active SMBH in the center of NGC 4561. 



The soft emission can be modeled as an absorbed thermal plasma with kT 
keV and non-solar abundances as follow: O = 0.33^^25 Na = 34t 8 , Si < 0.2, Fe 
and Ni < 0.4, with 1.0 corresponding to the Solar values. 



o sq +0 - 04 



0.12; 



-0.06 
-0.04 



No other model was capable of providing a good fit for the spectrum. It is quite in- 
teresting that non-solar abundances are needed for some elements. The nickel abundance is 
consistent with the solar value within the 3a contour. For oxygen, sodium, silicon and iron 
the abundances are not consistent with the solar values as can be seen in the contour plots 
(figuresEl [91 [10] and [11] for O, Na, Si and Fe respectively), they are sub-solar for oxygen, sili- 
con and iron but super-solar for sodium. It is important to note that when these parameters 
are kept frozen at solar value the model does not provide a good fit. Different abundances 
for different elements suggest that the observed thermal plasma is not well mixed in heavy 
elements or that the models we use are too simplistic to describe its physical conditions. 

Abundances for some elements in AGNs have been measured before, for exa mple in Mrk 



1044 and Mrk 279, wher e super-solar abundances fo r C, Ni, O and Fe were found ([Fields et al. 



2005. 


Fields et al. 


2007 





super-solar metallicities. Later studie s showed that quasars have super-solar abundances at 



all redshifts (IHamann fe Simonll2010l ) which would be consistent with the scenario of AGNs 
appearing after an important star formation event. Sub-solar abundance s, however, have 



been observed in the vicinity of AGNs like in the spectrum of NGC 1365 (jGuainazzi et al. 



20091 ). that shows sub-solar abundances of carbon, nitrogen, oxygen, neon, magnesium, sili- 
con and iron. Sub-solar abundance s were inferred al so for the sample of LLAGNs, LINERs 



and starburst galaxies presented in lPtak et al.lll999 



The temperature for the thermal plasma, kT = 0.591q'o5 keV, is similar to what has 
been found in similar objects like NGC 3367 and NGC 4 536 (kT = 0.64 ± 0.03 keV and 
kT = 0.58 ± 0.03 keV respectively, iMcAlpine et al.l 201 if ), and in agreement with what is 
expected for LLAGNs; kT « 0.4 - 0.8 keV (jHoll2008h ~ 



Levenson et al.ll200ll have analyzed ROSAT and ASCA data of Seyfert 2 galaxies with 
starbursts; they found that most of the soft emission, which is modeled with a thermal 
component, is produced by star formation. The median temperature of this component for 
their sample is about 0.6-0.7 keV. In their analysis, solar abundances were used because 
the quality of the data would not allow a measure of metallicity, but they do mention that 
low abundances are required in high-resolution spectra of starburst galaxies as shown by 
Dahlem et al.lll998l . Given the similarity of the temperature of the thermal component and 
the sub-solar abundances of NGC 4561 that we find, the soft X-ray emission we observe 
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could be produced by nuclear star formation. 



No variability was observed during the observation but there was a variation between 
the Chandra observation on 2006 March 15 and our XMM-Newton observation. From the 
count rate observed on Chandra (0.029 cts/s), PIMMS predicts a count rate of 0.072 cts/s for 
XMM PN using the thin filter and the whole PSF (~5 arcmin). We expect a more realistic 
value of about 78% of this rate, given the smaller (20" radius) PSF we used; it corresponds to 
0.056 cts/s. Our observed count rate of 0.282 ± 0.004 cts/s is 5 times higher than expected. 
This variation might have resulted from the change in the absorber column density. The 
variability also supports the AGN scenario. 

A lower limit on the black hole mass can be obtained assuming the black hole radiates 
at Eddington luminosity; for Lboi = 3 x 10 42 ergs s _1 , the BH must have a mass Mbh > 
2 x 10 4 M Q . This is an exciting discovery of a SMBH in a bulge-less galaxy. To our 
knowl edge, two bulge-less galaxies are known to host AG Ns from optical studi es: NGC 4395 



(Sdm; iFilippenko & Sargentlll989l IPeterson et all 120051 ) and NGC 1042 (Scd; IShields et al. 



20081 ). Both galaxies host a nuclear star cluster. IR spectroscopy with Sp itzer led to the 
discovery of AGNs in two more Sd galaxies, NGC 3621 rtSatvapal et aDl20071 ) and NGC 4178 
( ISatyapal et al.l 120091 ). These two also harbor a nuclear star cluster. Using XMM-N ewton , 



one m 



2 other AGNs in bulge-less galaxies have been confirmed by iMcAlpine et al.l 12011 
NGC 3367 (Sc) and the other in NGC 4536 (SABbc). We found nuclear X-r ay sources in 



M101 (type Scd), NGC 4713 (Sd), N GC 3184 (Scd), and NGC 4561 (Sdm) flGhosh et al. 



20081 . iGhoshl 120091 iGhosh et al.l 1201 if ) and argued that they are likely to be AGNs. Th e 



dwarf starburst galaxy Henize 2-10 is very likely to host an AGN too ( jReines et al.l 120111 ). 
In this paper we present conclusive evidence that NGC 4561 does in fact host an AGN. The 
discovery of SMBHs in the nuclei of these galaxies calls into question whether the masses of 
SMBHs are governed by bulge properties (§1). These results suggest that SMBH in bulge-less 
galaxies are far more common than what we previously thought, in clear disagreement with 
some models of BH growth. A key prediction of the models of Volonteri & Natarajan (2009) 
is that low-mass bulge-less galaxies today are unlikely to host nuclear black holes, which 
does not seem to be the case. While BHs grow through merger-driven processes, alternative 
tracks of BH growth must exist; secular process appear to play an important role. 



About 75% of late-type galaxies host nuclea r star-clusters (|Boker et al.l 120041 ) . Is the 
SMB H related to the m ass of the star-cluster then ( Seth et al. 20081 )? Or is it the dark matter 
halo ( jBaes et al.l 120031 )? Are the BHs in bulge- less galaxies the seed BHs at high redshift 
that did not grow? Answering these kind of questions is fundamental to our knowledge of 
BH-galaxy formation and co-evolution. 



6. Conclusions 



We present XMM-Newton spectrum of the Chandra-detected nuclear source in NGC 
4561 and show that it is an obscured AGN. The existence of nuclear SMBHs in bulge-less 
galaxies shows that BH masses are not governed by bulge properties. This calls into question 
several theoretical models (§1) of BH-galaxy co-evolution which are merger- driven; secular 
processes clearly play an important role in BH growth. 
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Fig. 1. — The EPIC-pn image of the source; the four panels show the different areas used to 
extract the spectrum when checking for possible contamination. To select the background 
spectrum the same 41" side square (in white) was used all the time. The different selection 
areas (in black) for the spectrum are: (1) 20" radius circle, (2) 26.65" radius circle, (3) 26.65" 
circle radius excluding a circle of the same size around the closest source and (4) a 26.65" 
radius semicircle opposite to the closest source. 
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Fig. 2. — The extracted spectrum and the model of an absorbed power law fitted to E>2.5 
keV. The bottom panel shows the contributions to \ 2 after extrapolating the model to the 
entire energy range. Significant absorption is needed to fit the hard end of the spectrum 
properly. 
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Fig. 3. — The data fitted with an absorbed power-law plus the diffuse thermal plasma model 
( "mekal" ) with solar abundances and a black body. Contributions to x 2 are presented in the 
bottom panel and show that this model does not fit the data well. 
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Fig. 4. — Same as Fig. 3, but with a constant non-solar abundance (each component is 
shown with a dotted line and the solid line is the sum of all components). This model does 
not fit the data well either. 
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Fig. 5. — Same as Fig. 4, but with variable abundances. This is the best-fit model. 
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Fig. 6. — The spectrum fitted with an absorbed power-law plus an ionized reflection compo- 
nent. Contributions to x 2 are presented in the bottom panel. This model does not provide 
a good fit. 
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Current Theoretical Model 
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araya S-Deo2011 17:49 

Fig. 7. — Theoretical model of the best fit spectrum showing the different components 
(absorbed power law, absorbed thermal plasma and black body). 

Confidence contours: Chi-Squared 
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araya 13-Dec-2011 16:07 

Fig. 8. — x 2 contours for the abundance of oxygen and temperature of the plasma in the 
best fit model (Solar abundance is 1). Sub-solar oxygen abundance is clearly required. 
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Confidence contours: Chi-Squared 
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Fig. 9. — x 2 contours for the abundance of sodium and temperature of the plasma in the 
best fit model (Solar abundance is 1). Sodium abundance appears to be super-solar. 
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Fig. 10. — x 2 contours for the abundance of silicon and temperature of the plasma in the 
best fit model (Solar abundance is 1). Sub-solar silicon is clearly indicated. 
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Confidence contours: Chi-Squared 
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Fig. 11. — x 2 contours for the abundance of iron and temperature of the plasma in the best 
fit model (Solar abundance is 1). Again, sub-solar iron abundance is required by the fit. 
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Table I. Spectral Models 



Model ID 



Model 



Component 



Parameters* 



Xl (dof) 



absorbed power-law and thermal power-law 

plasma with solar abundance mckal 

absorbed power-law and thermal power-law 

plasma with non-solar abundance mekal 

absorbed power-law and power-law 

reflection rcflionx 

absorbed power-law power-law 

and vmckal with solar abundance vmekal 

absorbed power-law and power-law 

vmckal with non-solar abundance vmekal 



N H = 1.02 ± 0.05 
N H = 0.86 ± 0.05, kT = 0.11 ± 0.01 

Nf/ = 1.06 ± 0.05 
N H = 0.82 ± 0.05 , kT = 0.12 ± 0.01, abund= 0.3 ± 0.3 



: 0.23 



+0.03 



Ntf = 1.6 ± 0.2 
Fc/solar = 4.5+2:5, X * 



1300 



+ 700 



N H = 1.39 ±0.06 
N H = 0.0+ 03 , kT = 0.70 ± 0.01 

N H = 1.9l u2 
JV„=0.07l°;°§,fcT = 0.59±°;gS 
O = 0.33±g 12 % Na = 34+^ 
Si = 0±°"^ Eb = b.l2±°;<» Ni = 0±°"' 



2.35 (112) 
2.27 (111) 
1.39 (113) 
2.35 (112) 
1.11 (107) 



*In all the models: Galactic absorption was added with N^f= 2.11 X 10 20 cm 2 and the photon index of the power-law is T = 
2.46 (obtained fitting only for E > 2.5 kcV). These two parameters were always frozen when fitting. All the models also have a 
blackbody with kT=0.2 kcV and flux F(0.2-10 kcV)=2.5xlO~ 14 accounting for the source B. 

**A1 column densities are in units of 10 22 cm" 2 , temperatures in keV. 



